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Third, a discrete symmetry group and its fixed points which are Maxwell points of of timelike 
and spacelike normal extremals, are described. An estimate for the cut time (the time of loss 
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1 Introduction 

A sub-Riemannian structure on a manifold M is given by a smoothly varying distribu¬ 
tion D on M and a smoothly varying positive definite metric g on the distribution. The 
triple (M, D,g) is called a sub-Riemannian manifold, which has been applied in control theory, 
quantum physics, C-R geometry and other areas. Some efforts have been made to generalize 
sub-Riemannian manifold. One of them leads to the following question: what kind of geo¬ 
metrical features the mentioned triple will have if we change the positive definite metric to an 
♦Supported by NSFC (No.ll071119, No.ll401531); NSFC-RFBR (No. 11311120055). 
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indefinite nondegenerate metric? It is natural to start with the Lorentzian metric of index 1. In 
this case the triple: manifold, distribution and Lorentzian metric on the distribution is called 
a sub-Lorentzian manifold by analogy to a Lorentzian manifold. For details concerning the 
sub-Lorentzian geometry, the reader is referred to HZ]. To our knowledge, there are only few 
works devoted to this subject (see HUlIZlIISlIISlEniES]). In HU, Chang, Markina, and Vasiliev 
systematically studied geodesics in an anti-de Sitter space with a sub-Lorentzian metric and 
a sub-Riemannian metric respectively. In [19], Grochowski computed reachable sets starting 
from a point in the Heisenberg sub-Lorentzian manifold on M^. It was shown in |25] that the 
Heisenberg group H with a Lorentzian metric on possesses the uniqueness of Hamiltonian 
geodesics of time-like or space-like type. 

The Engel group was first named by Cartan [9] in 1901. It is a prolongation of a three 
dimensional contact manifold, and is a Goursat manifold. In laisizi A. Ardentov and Y. 
Sachkov computed minimizers on the sub-Riemannian Engel group. In the present article, we 
study the Engel group furnished with a sub-Lorentzian metric. This is an interesting example of 
sub-Lorentzian manifolds, because the Engel group is the simplest sub-Lorentzian manifold with 
nontrivial abnormal extremal trajectories, and the vector distribution of the Engel group is not 
2—generating, its growth vector is (2,3,4). We first prove that timelike normal extremals are 
locally maxiziming. Second, we study the normal extremals. By using the Pontryagin Maximum 
Principle, the subsystem for costate variables of normal Hamiltonian system is reduced to 
equations similar to pendulum. Expressions of timelike, spacelike and lightlike normal extremal 
trajectories are obtained in terms fo Jacobi functions. Third, a discrete symmetry group and 
its fixed points, which are Maxwell points, are described. An estimate for the cut time (the 
time of loss of optimality) on extremal trajectories is derived on this basis. 

The structure of this paper is as follows. Section 2 contains some preliminaries as well 
as definitions of sub-Lorentzian manifolds and the Engel group. In Section 3, we prove that 
timelike normal extremals are locally maxiziming. In Section 4 we get a description of light- 
like normal extremal trajectories. In Section 5, we find a parametrization of timelike normal 
extremal trajectories and describe the exponential map providing the parametrization for all 
the timelike normal extremal trajectories, then we describe symmetries of the exponential map 
and investigate the corresponding Maxwell points, which are fixed points of these symmetries. 
In Section 6, a parametrization of spacelike normal extremal trajectories are obtained, and the 
exponential map providing parametrization for all the spacelike normal extremal trajectories is 
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described, then we describe symmetries of the exponential map and investigate the correspond¬ 
ing Maxwell points, which are fixed points of these symmetries. On the basis of this study, we 
prove upper bounds on the cut time. 

2 Preliminaries 

A sub-Lorentzian manifold is a triple (M, D,g), where M is a smooth n—dimensional man¬ 
ifold, ID is a smooth distribution on M and o' is a smoothly varying Lorentzian metric on D. 
For each point p G M, a vector v G Dp is said to be horizontal. An absolutely continuous curve 
7 (t) is said to be horizontal if its derivative 7 ^(t) exists almost everywhere and lies in D^^^y 
A vector v G Dp is said to be time-like if g{v, v) < 0; space-like if g{v, ?;) > 0 or r; = 0; 
null (lightlike) if g{v^ u) = 0 and v ^ 0] and non-space-like if g{v, v) < 0. A curve is said to 
be time-like if its tangent vector is time-like a.e.; similarly, space-like, null, and non-space-like 
curves can be defined. 

By a time orientation of (M, D, g), we mean a continuous time-like vector field on M. From 
now on, we assume that {M,D,g) is time-oriented. If A is a time orientation on {M,D,g), 
then a non-space-like v G Dp is said to be future directed if g(v,X{p)) < 0, and past directed if 
g{v,X{p)) > 0. Throughout this paper, “f.d." stands for “future directed", “t." for “time-like", 
and “nspc." for “non-space-like". 

Let v,w € D he two nspc. vectors, we have a following reverse Schwarz-inequality |24| : 

|5r(r;, w)| > Hull • ||u;||, (1) 

where ||u|| = ^/\g(vyv)\, the equality holds if and only if v and w are linearly dependent. If v is 
a nonspacelike f.d. vector and w is a nonspacelike p.d. vector, it is easy to see that g(v, w) > 0. 
We introduce the space of horizontal nspc. curves: 

^7(t) = {7 : [0, 1 ] -t M| 5 r( 7 '(f), 7 '(t)) < 0,7'(f) G D^y^ for almost all t G [0,1]}. 

The sub-Lorentzian length of a horizontal nspc. curve 'y{t) is defined as follows: 

^( 7 ) = [ \\'y\t)\\dt, 

Jo 

where ||7^(t)|| = is computed using the Lorentzian metric on the horizontal 

spaces D^yy We use the length to define the sub-Lorentzian distance du(qi,q 2 ) with respect 
to a set 1/ C M between two points qi,q2,& U- 
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du{qi,q2) 


( sup{/(7),7 G ^[ 7 ( 51 , 92 )} if Hu{quq2)^^ 

0 if Hu{qi,q2) = 0, 

where Hu[qi,q 2 ) is the set of all nspc.f.d curves contained in U and joining qi and ^ 2 - 

A distribution D C TM is called bracket generating if any local frame {Aij}i<j<r for D, 
together with all of its iterated Lie brackets \Xi,Xj\, [Xi, [Xj, X^]], • • • span the tangent bundle 
TM. Bracket generating distributions are also called completely nonholonomic or distributions 
satisfying Hdrmander’s condition. 

Theorem 1. (Rashevsky-Chow) Fix a point q G M. If a distribution D C TM is bracket 
generating then the set of points that can be connected to q by a horizontal curve is the component 
of M containing q. 

By the Rashevsky-Chow Theorem, we know that if D is bracket generating and M is 
connected, then any two points of M can be joined by a horizontal curve. 

In the Lorentzian geometry the lengh minimizer problem is trivial because arbitrary two 
points can be connected by a piecewise lightlike curve of zero length. However, there do exist 
timelike curves with maximal length which are timelike geodesics |24) . Thus it is natural to 
study sub-Lorentzian length maximizers. 

A nspc. curve is said to be a maximizer if it realizes the distance between its endpoints. 
We also use the name [/-extremal for a curve in U whose each suitably short sub-arc is a 
[/-maximizer. A nspc. curve will be called geodesic if it is a local maximizer. 

The Hamilton function H : T*M ^ M he defined as 

V 

= <A,Xi(x)>2+1^<A,W(x)>2, {x,\)eT*M. (2) 

^ i=2 

Let U C M be an open neighborhood, ip : U cM—J-Risa smooth function. 

Definition 2.1. The horizontal gradient XhT of p is a smooth horizontal vector field on U 
such that for each p £ U and v £ H we have dyp{p) = g(V hp{p),v). 

Locally we can write 

r 

XhP = -{dx,p)X^ + Y.{dxMX., (3) 

i=2 

Definition 2.2. A normal extremal trajectory in the sub-Lorentzian manifold (M,D,g) is a 
curve 7 ([) : [a, 6] —)■ M that admits a lift F : [a, b] —>■ T*M which is a solution of the Hamiltonian 
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system with the sub-Lorentzian Hamiltonian H{x,X). In this ease, we say that T{t) is a normal 
lift ofj{t). 


Now we introduce the Engel group E with coordinates q = {xi,X 2 ,y, z) G The group 
law is denoted by © and is defined as follows: 


{xi,X2,y,z) 0 {x[,X2,y',z') 

= Tx! + Xi,X2 + X2,y + y + ^1^2 ^ + X2) + xiy' + + x[) ] . 


A vector field X on El is said to be left-invariant if it satisfies dLqX{e) = X{q), where Lq 
denotes the left translation p ^ Lq{p) = qQp and e is the identity of E. This definition implies 
that any left-invariant vector field on El is a linear combination with constant coefficients of the 
following vector fields: 




^3 


d X2 d d xi d x1 + X 2 d 

dxi “I dy' ^ dx 2 2 2 dz' 

d d d 

oy oz oz 


(4) 


The distribution D = span{Xi,X 2 } on E satisfies the bracket generating condition, since 
X 3 = [Xi,X 2 ], X 4 = [XijXs], The Engel group is a nilpotent Lie group, since [Xi,X 4 ] = 
[X 2 ,X 3 ] = [X 2 ,X 4 ] = 0. We define a smooth Lorentzian metric on D by 


ff(Xi,Xi) = -l, g{X2,X2) = l, g{XuX2) = Q. (5) 


3 Local maximizing property of normal extremals 

In this secton, we prove that timelike normal extremals are locally maximizing. 

Theorem 2. Let ^{t) be a t.f.d. normal extremal defined on an interval [a,b]. Then for eaeh 
t G [a,b], there exists a neighborhood U of ^{t) sueh that the restrietion of'j to U is a unique 
maximizer eonneeting its endpoints. 

Proof Let T{t) = {'y{t), X{t)) G T*M be a normal lift of 7 (t), then we can assume that 
H{'y,X) = — G {a,b). To show that 7 is locally maximizing, we need only to prove 

that for every c G {a,b), there exist J = (ti,t 2 ) © (O; ^ small interval of c, such that the 
restriction of 7 (t) to J is maximizer. Let pc = 7 (c) and Ac = A(c). Pick a smooth hypersurface 
S' C M of dimension n — 1 in M passing through pc such that Ac vanishes on Tp^S. Let A be 
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a smooth one-form on an open neighborhood of Pc such that X{pc) = A and A(p) annihilates 
TpS and H{p,X) = — ^ for all p G 5 Pi 17. 

Let Lp = ( 7 p,Ap) be the solution of f(t) = ^(r(f)), r(c) = (p, A), p G 5 H n. It is clear 
that r = Lp^. Since 7 ^ Tp^S, by the Implicit Function Theorem, we can get, if e > 0 is small 
enough and IF is a sufficiently small open neighborhood of Pc in S, that there exist a map 

u :{c — e,c + e) X W ^ M, 
it,p) lp{t), 

that maps (c — e, c + e) x IF diffeomorphically onto an open neighborhood U of pc in M. Define 
a smooth function (p : U ^ R and a 1-form 0 on 1/ by letting ‘p{x) = t and 6{x) = \p{t) if 

X = -fp{t). 

We will prove that dp = —6. Let X be a vector field on U such that 

x( 7 (t)) = d7r,(rp(t))t(rp(t)), 

then we know that jp{t) = X{'yp{t)). Let Hx he the Hamiltonian lift of W, we can show that 
fx(rp(t)) = t(rp(t)), VpGiF, te{c-E,c + e). 

Indeed, let ai = {9{x),Xi{x)) for i = 1, ..., /c, if x G U, then 

k 

X = —aiXi + aiXj. 

i=2 

Since H is constant along integral curves of then 

H{hp{t)) = V(t,p) G (c-e,c + e) X IF, 

so we have that 

k k 

— ai{x) + '^^af{x) =—l, —ai(x)ai(x)ai(x)^(x) = 0, Vx G C/. ( 6 ) 

i=2 i=2 

Since 

V ^ V 

t(x, A) = -(A, Xi):^(x, A) + J^(A, Xi)H^Xx, A), (7) 

i=2 

in particular 

k 

ilpit), Xp{t)) = -ai{jp{t))I^Xiilpit), Xp{t)) + ^ai{jp{t))I^Xihpit)^ Xp{t)). 

i=2 
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On the other hand 


'^x{x,X) = -ai{x)'^xi{x, X) + y^^ai{x)'^x,{x,X) - {X,Xi)at{x) + (A,Xj)^(x). 

i=2 


Therefore 


_^ _ K 

Hxilpi Xp) = —ai{'^p)^Xi{lp-,Xp) + Xii'Jpi Xp) 

i=2 

k 

— {Xp, Xi('jp))ai{'yp) + 'y ^^ {Xp, Xi{^p)) ai{^p 

i=2 

k 

— ®i(7p)^Xi i'lpi Xp) + ^ Xjilpi Xp') 

i=2 

k 

- ai(7p)^(7p) + ^ai{jp)'(^i{jp) 

i=2 

k 

— ®i(7p)^Xi {'Ipi Xp) + ^ Xjilpi Xp') 

= ^(7p, Ap). 


i=2 


So, for p G tT, Tp is also the solution of 


t{t) = tx{T{t)), T{c) = {p,X{p))- 


( 8 ) 


Let 4>t be the flow oi X onU for t G (c — e, c + e) and let (t)c{p) = P for p G U. For p G W, 
we have 7 p(t) = 

For each t G (c — e, c + e), let IFt = {x G C/ : p{x) = t}. Then Wt is a smooth hypersurface, 
and Wt = 4>t{W) and T^Wt = d(t>t{p)TpW for x = 7 p(t). 

We now show that dp = —6 on U. Let x G U and p G W,t G (c — e, c + e), be such that 
X = 'Jpit), then X G IFt. Let w G T^M and v G TpM be such than w = d(j)t{p)v. Define 
h{s) = d(j)s{p)v, then h{s) G T^^(^g-jWs and h{t) = w. It is easy to see that h{s) satisfies the 
following equation: 

f) X 

h{s) = 

Since Ap satisfies the Hamiltonian equations, so 

7(») = -^(>(>)) = 

hence the fucntion s i—)• (Ap(s), h{s)) is constant on (c — e, c + e). 
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If tj; G T^Wt, then v G TpW. Since fix) is constant on Wt, then {d(p{x),w) = 0. Since 
{Xp{c),v) = 0, we have that 

{9{x),w) = {Xp{t),h{t)) = (Ap(c), h{c)) = (Ap(c), z;) = 0 = -{dip{x),w). (9) 

If u) = 7 p(t), then v = 'jpic). Since {dip{x),'yp{t)) = 1, we have that 

= {Xp{t),h{t)) = {Xp{c),h{c)) = -1 = -{d(p{x),w). (10) 

Since T^M = TxWt © { 7 p(t)}, we conclude that 9{x) = —d(p{x). 

By ([^, (|^ and ( p©] ), we can get that 

( 11 ) 

Since 'y{t) is a t.f.d curve, V is a p.d. vector field on U. 

Since H{x,9{x)) = —^ and ||0(x)|| = —1, we have shown that tp satisfies the Hamiltonian 
equation H{x,dip{x)) = —^ and ||(i(^(x)|| = —1. 

Next, choose ti and t 2 in the interval (c — e, c + e). Now, let 7 : [0, a] be a t.f.d. curve in U 
with 7 ( 0 ) = 7 (ti) and r/(a) = 7 (t 2 ); then we have 

^(7l[7.t2]) =t2-ti= <^( 7 (^ 2 )) - 

pOL pCt 

= p{r]{0)) - ip{T]{a)) = {dip{r]{s)))ds= g{f],VhP) 

Jo Jo 

> [ ||?7(s)||ds = L(7 |[o,„]). (12) 

Jo 

By the reverse Schwarz inequality, L{'y) = L{r]) holds if and only if rj can be reparameterized 
as a trajectory of —’VhP- This completes the proof. □ 

Remark 1. Let ^{t) be a t.p.d. normal extremal. We also can prove that ^{t) is a the unique 
local maximizer by the above method. 


4 Sub-Lorentzian extremal trajectories 


In this section, we calculate sub-Lorentzian extremal trajectories by applying the Pontryagin 
maximum principle. We can assume that the initial point is the origin by invariance under left 
translations of the Engel group, i.e., xi(0) = X2(0) = y(0} = z(0) = 0. Let us introduce the 
vector of costate variables ^ = (^O) ? 2 , ^ 3 , ^ 4 ) and define the Hamiltonian function 

X 2 U 1 


H{i,q{t),u) = go +gwi +g2^^2 


^2 I /ya2 

+ U^^V^U2. (13) 


2 


2 


2 





From the Pontryagin maximum principle for this Hamiltonian function we obtain a Hamiltonian 
system for the costate variables: 


6 = -H,, 


^ 3^2 

2 


UxiU2, 


i = -h.2 


2 


CiX2U2, .^3 


^4 = 0 , 


(14) 


the maximality condition 


maxH{^{t),q{t),u{t)) = H{^{t),q{t),u{t)), < 0, 


(15) 


where u{t),q{t) is the optimal process, and the condition ^{t) 7 ^ 0 of nontriviality of the costate 
variables. 

For sub-Lorentzian extremal trajectories, there are abnormal extremal trajectories (i.e. .^0 = 
0) and normal extremal trajectories (i.e. ^0 = ~1)- Abnormal extremal trajectories were 
described in HHQ] 


Now we look at the normal case ^0 = It follows from the maximality condition (15) 


that = 0. Hence 


r. X2^3. . , , Uixj + x^) 

Ul = -i^l - U2=^2 + — + - - -. 


(16) 


Let hi = {^,Xi),i = 1, • • • ,4, be the Hamiltonians corresponding to the basis vector fields 
Ai,X 2 , A 3 , A 4 in the tangent space TqE and linear on the fibres of the cotangent space T*E: 

2 2 

hi = Cl —^^3) h2 = C 2 + H— ^ ^4, h3 = C 3 + xiCi, hi = ^ 4 . (17) 


So ui = —hi and U 2 = h 2 - 

According to the Pontryagin maximum principle, we get the following Hamilton equations 
in the new coordinates: 
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( 18 ) 


"1 = % = -<fi - 


dH xi^ Xi + Xn^ , 

^2 = =6+ 2 ^ 3 + 2 ^4 = /i2, 

(19) 

y = ^ = + ^2y = ^{xih2 + X2hi), 

(20) 

dH X? -|- Xn , 

2 

(21) 

^1=6“ y^3 = -h2h3, 

(22) 

X 2 

h 2 =f .2 + y 6 + Uixixi + X 2 X 2 ) = -hihs, 

(23) 

1 

II 

•H 

II 

CO 

(24) 

II 

0 

(25) 


Associated with the expression of H, we conclude that a sub-Lorentzian extremal is timelike 
if < 0; spacelike if H > 0] lightlike if H = 0. 

For the lightlike case, by the definition, we have H = \{—h\ + /i|) = 0, thus /12 = ±/ii. If 
h 2 = hi, then lightlike trajectories satisfy the ODE: 

7 = -hi(Xi-X2), 


that is, they are reparameterizations of the one-parametric subgroup of the field Xi — X 2 . We 
assume 7 = Ai — A 2 , so 

xi = 1, X2 =-1, y = -^{xi + X2), z =-^(xj + xl), 

thus 


1 s 

Xi=t, x 2 = -t, y = 0, z = --r. 
If /i 2 = —hi, similarly, we obtain 

1 o 

xi=t, X 2 = t, y = 0, z = -t . 
In conclusion, we get the following theorem: 


Theorem 3. Lightlike horizontal extremal trajeetories starting from origin are reparameteriza¬ 
tions of the eurves: 

1 o 

xi = t, X 2 = ±t, y = 0, z = ±-t . 

O 

Timelike and spacelike normal extremal trajectories are considered in the following two 
sections to discuss them. 


10 



5 Timelike normal extremal trajectories 


In the timelike case H = + h"^) < 0, we consider extremals on the level surface 

{H = —1/2} and introduce coordinates (0,c,a) on this surface as follows: 

hi = ±cosh 0 , /i 2 = sinh 0 , /13 = c, /i 4 = a. 


Notice that we have the following symmetry of the Hamiltonian system (18)-(25): 


: (hi, h 2 , / 13 , hi, xi,X 2 ,y, z) ^ (-hi, h 2 , -h^, hi, -xi,X 2 , -y, z). 


(26) 


So we consider case hi = cosh0 > 0 without loss of generality in the sequel. 

In the variables (0,c,a,xi,X2,y,z) on the level surface {H = —1/2} the Hamiltonian sys¬ 


tem (18)-(25) assumes the following form: 


xi = — cosh 0 . 

(27) 

i ;2 = sinh 0 . 

(28) 

X 2 cosh 9 + xi sinh 9 

y= 2 

(29) 

X 1 +X 2 

2 

(30) 

9 = —c. 

(31) 

c = —a cosh 0 , 

(32) 

d = 0 . 

(33) 


Note that the subsystem for the costate variables reduces to the equations 

6 = a cosh 0 , d = 0 , 


(34) 


whose phase portrait for a = 1 and a = — 1 is given in Fig. 

Let us introduce the energy integral: 

h2 2 

E = — / 12/14 = -^— asinh^, E = / 13/13 — / 14/12 = 0. (35) 

The set of all normal extremal paths is parametrized by points of the set 


C = T*^M n {H = — 1 / 2 , hi > 0} = {(hi, h 2 , h^, hi) G | h{ — h^ = 1, hi > 0} 
= {(h,c,a)GM='}. 
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Figure 1: Phase portrait of vertical subsystem (34) for C~^ and C (timelike case) 


The set C has the following decomposition into invariant sets of equations (31)-(33): 


C = CSUC°UC+UC-, X = {9,c,a), 
Cq = {A G C I c = 0, a = 0}, 

C° = {AgC| c^0,a = 0}, 

C+ = {A G C I a > 0}, 

C- = {A G C I a < 0}. 


Consider the general case a 7 ^ 0. Trajectories of vertical subsystem (31)-(33) for and 
C~ are symmetric (see Fig. [^. This symmetry (see Section 5.2) has the following description: 

: (a, c, 6, xi,X 2 ,y, z) 1 -^ {-a, -c, -0, xi, -X 2 , -y, -z). (36) 


Thus, it is enough to integrate the Hamiltonian system in the case due to the sym¬ 
metry e^. We introduce new coordinates {ip,E,a) in the subset which rectify the vertical 


subsystem (31)-(33): 


as = 






c = —2aesc(8e(p) dn(8e(p), 




sinh 0 = 


2as^^l — k‘^[l + 


K K\ 

’ / ■ 
« ae/ 


cn 


a cn2(as(p) 


cosh 0 = 


2m^( 1 — /c^(l — 


cn 


a cv?{?e^p) 
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where sn^, cn'0, dn are elliptic Jacobi functions with modulus k, and sc'ip = 


snip 

cmp' 


K(fc) = 


dt 


is the complete elliptic integral of the hrst kind. 


/o s/l — sin^ t 

Immediate differentiation shows that in these coordinates the subsystem for the costate 


variables (31)-(33) takes the following form: 


(^ = 1 , E = 0, d = 0 , 


so that it has the solutions 


(p{t) = (ft = 'fo + t, E = const, a = const. 


(37) 


(38) 


5.1 Exponential mapping for timelike normal extremals 

Denote arguments of Jacobi functions pjo = G (~K)K) to describe the 

exponential mapping in general case a ^ 0: 


2 se 

xi{t) = ■^{sc'ipodnijjQ - sciptdnipt), 


\a\ 


X2{t) = 


4 ae(^as(l - k'^)t - (£(?/)*) - E(V’o))) - |aki(t) 


a 


y{t) =-r-rf fe^(cn^ V't “ cn^ V’o) + (1 - k‘^){nc^'ipt - nc^ 

a\a\ \ 

* ^dn'0osc'0o+ dn'i/’tscV't^X2(t), 


(39) 

(40) 

(41) 


(X 2 (t))^ 486^ 


2 ./dn sn V't dnV'osn'^o 


+ fe^(cn?/)idn?/)t snV'f - cnipo dnipo sm/jo) - 2{E{'ijjt) - E{ipo)){2k‘^ - 1) 

28edn'!/’osnV’o(^^^^-]^]^^a;2(t) - lxi{t)x2{t) - 


cn ^po\a\ 


2ce^{2k‘^ — l)xi(t) 
3 q;|q:| 


(42) 


rip 


where E(^/;) = / dn^ t dt. 

Jo 

Projections of extremals to the plane (xi, X 2 ) with \a\ = 1 and ■00 = 0 are shown in Fig. 
left. Notice that since 'ipt G (~ K, K) we have an upper bound for the time parameter 

K{k) - 00 


t ^ tsupr(^) — 


ae 


(43) 


When t —7- tsupr(A), solutions (39)-(42) explode, i.e., tend to inhnity. 
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Figure 2: Projections of extremals to the plane {xi,X 2 ) for C^, and Cq (timelike 
case) 


For the subset we have 


a = 0, E = 


c = const, 6 = 6o — ct, 


Xl{t) = 

X2{t) = 

y{t) = 


smh( 0 o — ct) — sinh^o 

? 

c 

cosh 00 — cosh( 0 o — ct) 
c 

sinh(ct) — c t 


2c2 


z{t) = - 


4sinfo (^) sinh (30o — — 3(sinh(ct) — ct) sinh0o 


6c^ 


Projections of extremals to the plane {x\,X 2 ) with |c| = 1 (hyperbolas) are shown in Fig. 
center. 


If A G Cq, then 


q; = c = £^ = 0, 0 = const, sinh 0 = sq, cosh 0 = cq, 

Xi{t) = -Co t, 

X2{t) = So t, 
y{t) = 0, 

z{t) = t\ 

Projections of extremals to the plane {xi,X 2 ) are straight lines (see Fig. right). For A G 
U Cq we have no upper bound for the time parameter: tsupr(A) = +oo. 
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So we obtained a parameterization of the exponential mapping 


Exp : iV ^ M = 


Exp(A,t) = q{t), 


N — {(A,t) £ C X I t £ (O, tsupr(A))}, 


(44) 

(45) 


in the timelike case. It maps a pair (A, t) to the point of the corresponding extremal trajectory 
q{t). 

Further we investigate discrete symmetries of the exponential mapping and on this basis 
find estimates for the cut time on extremals. 


5.2 Discrete symmetries of exponential mapping 


Subsystem (31 )-(33) for costate variables of the normal Hamiltonian system has symmetries 


which preserve the direction field of this system. We already described one of them (36) 


Let us define the action of symmetries e® on the set of trajectories of the vertical subsystem 
with preservation of time direction. Denote a smooth curve 

7 = [{0{t),c{t),a) I t £ [0,r]| C C. 

Define the action of symmetries on these curves (see. Fig. [^: 

: 7 ^ 71 = (t), (t), a^) \ t £ [0,r]| = {( - ^(t),-c(t),-a)}, 

^2 = {{0\t),c^it),a^) \t£[0,T]] = {{9iT-t),-c{T-t), a)], 

: 7 ^ 73 = | (<9^(i), c^(i), I ^ e [0, r]| = | ( - e{T - t),c{T - t), -a)|. 


It is obvious that if 7 is a solution to the vertical subsystem (31)~(33), then 77 i = 1,2,3, is a 
solution as well. 

Now consider the group of symmetries G = {Id,o = Z 2 x Z 2 . The 
symmetry preserves the direction of time, while and change it. 

Continue the action of symmetries from the vertical subsystem to solutions of the Hamilto¬ 
nian system of PMP 


9{t) = —c{t), c{t) = —acosh 9 (t), d = 0, 
q{t) = - cosh 9{t)XI {q{t)) + sinh 9{t)X 2 {q{t)), 


(46) 

(47) 


in the following way: 


z{t),a,q{t)) I t G [0,r]| |( 6 '®(t),c®(t),a®,g®(t))|f G [0,T]|, (48) 
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Figure 3: Action of symmetries on trajectories of vertical subsystem (timelike case) 
where q{t) = z{t)), t £ [0,T], is a geodesic and 

q^t) = {x\{t),xl{t),y\t),z\t)), t £ [0,r], i = 1,2,3, 
are its images under the action of the symmetries e*. 

Lemma 1 . The symmetries e* map trajeetories on the plane (xi,X2) in the following way: 


II 

xl{t) = -X 2 {t), 

xl{t) = xi{T) - xi{T - t), 

X^it) = X 2 (T) — X 2 {T — t) 

xl{t) = xi{T) - xi{T - t), 

xl{t) = X2{T -t) - X2{T) 


Proof. Via direct integration, for instance, for we verify that 

^i(^) = [ { —cosh9{T — s))ds = f [ —cosh 9{r))dr = xi{T) — xi{T — t), 

Jo Jt-i 

x\{t) = f sinh 6{T — s)ds = f sinh 9{r)dr = X 2 {T) — X 2 (T — t). 

Jo Jr-t 

□ 

Lemma 2. The symmetries e* map endpoints of geodesies q = {xi, X 2 ,y, z) to endpoints of 


16 



= {x\,X 2 ,y'', in the following way: 


xl(r) = xi(r), 

x\{T) = -X2{T), 

y\T) = -y{T), 

z\T) = 

-ziT), 

x?(r) = xi(r), 

xUT)=X 2 {T), 

y\T) = -y{T), 

z\T) = 

ziT)-xi{T)y{T), 

x?(r) = xi(T), 

xl{T) = -X2iT), 

y\T) = y{T), 

isi 

CO 

3 

II 

-z{T)+xi{T)yiT) 


Proof. Lemma gives us the expressions for x\{T) and x\{T). The expressions for other 
variables are obtained by integration. For example, for we have 

y^T) = £ + cost>(T - t)dt 

+ ?iT) Tsinhon-t)*- r d, 

2 7 o Jo 2 


X2(r)xi(r) xi(r)x2(r) 


-y(T) = -y{T), 




(i' 2 (r))° {x,{T)f - 2xi(T)xi{T -t) + (xi(T - t))" 


(x 2 (T})\ (x,(T)fx 2 (T) /-T Mr))' 

+ ^- xi{T) / Xi(s)x2(s)ds + 2;(T)- 


6 


6 


- ^i(r) f2/(r) + ) + HT) = z(T) - xi(T)y(T). 


□ 


We define the action of e* in the preimage N of the exponential mapping by restricting the 
action to initial point of trajectory of vertical subsystem: 


e^:C^C, e\e,c,a) = {e\c\a^), (49) 

{6^,c^,a^) = {—6,—c,—a), (50) 

(0^, c^, a^) = {6, —c, a), (51) 

{9^,c^,a^) = {-9,c,-a), (52) 


in the following way: 

e^(A,t) = (e^(A),t), 

e*(A, f) = (e* o e^^”(A),t), i = 2,3, 

where H^} = —cjg — acosh^J^; G Vec(C) is the vertical part of Hamiltonian vector field. 
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We define the action of e* in the image of exponential mapping M by restricting the action 


to endpoints of extremals (see 

Lemma : 


:M ^ M, 

£\q) = £\xi,X2,y,z) = g® = {x\,xl,y\z^), 

(53) 


= (xi, -X2, -y, -z), 

(54) 

{xl,xl,y‘^,z‘^) 

= (xi, X2, -y, z - Xiy), 

(55) 

{xl,xl,y^,z^) 

= (a:i, -X2, y, xiy - z). 

(56) 


Since the action of e* in the domain N and the image M of the exponential map are induced 


by the actions of symmetries (48) on trajectories of the Hamiltonian system (46)-(47|, we have 
the following result. 


Proposition 1. The mappings = 1,2,3, are symmetries of the exponential mapping, i.e., 


(e® o Exp)(0, c, a, t) = (Exp o e®)(0, c, a, t), 


{6, c, a, t) G N. 


5.3 Maxwell points of timelike normal extremals 

A point q{t) on a geodesic is called a Maxwell point if there exists another geodesic q{s) ^ 
q{s) for which q(t) = q{t),t > 0. It is known that a geodesic cannot be optimal after a Maxwell 
point. In this section we compute the Maxwell points corresponding to some symmetries e®. 
On this basis we derive estimates for the cut time along geodesics 

Aut(A) = sup{t > 0 I Exp(A, s) is optimal for s G [0,t]}. 

We define Maxwell sets in the preimage of Exp corresponding to the symmetries e®: 

MAX® = {(A,t) G I AV A, Exp(A^^) = Exp(A,t)}, (57) 

X = {9,c,a), A® = e®(A). 

It follows from Proposition that the equality Exp(A®,t) = Exp(A,t) is equivalent to 
^*(^(0) = Therefore we get a description of fixed points of the symmetries e® in the 

image of the exponential mapping. 

Lemma 3. 

1. e^{q) = q X 2 = 0, y = 0, z = 0. 
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2. e^{q) = q 

3. e\q) = q 


y = 0. 


X2 = 0, z = —. 


Proof. Follows from definition (53)-(56) of action of the symmetries in E. 


□ 


Now we compnte fixed points of symmetries in the preimage of exponential mapping, which 
is necessary for description of Maxwell sets. Introdnce the following coordinates in the sets 
Ar± = {(A,t) e X M+ I t G (0,4upr(A))}: 


r = 


V’t + V’o 


P = 


- -00 ast 


2 ’ ^ 2 2 
The parameter r corresponds to the middle point of a trajectory. 


(58) 


Lemma 4. 

1. = X ^ 

2. X'^ = X ^ 

3. X^ = X ^ A E U 0 + = 0 

Proof. Follows from the definition of action of reflections in the domain of exponential map¬ 


AeC°, 3 = 0. 

XgC^,t = 0 or A ECO. 


ping (49)"(52). 


□ 


Theorem 4. Let (A,T) E N and q{T) = Exp(A,T). Then 

1. (A, T) E MAX^ ^ X 2 {T) = y{T) = z{T) = 0 and if X £ then 9^0. 

(A,r)EMAX2 ^ AEC±,y(r) = 0,r/0 or X £ C°,y{T) = 0. 

3. (A, T) £ MAX^ ^ 

Proof. Follows from definition (57) of Maxwell sets, Lemmas and 


X 2 = 0, z = and if X £ U Cq, then 9 ^ . 


ct 


□ 


The eqnations e* (q) = q define snb-manifolds in with dimension from 3 to 1, exponential 
mapping translates the corresponding Maxwell sets into these snb-manifolds: 


Exp(MAX*) c{q£M \ e\q) = q}, i = 1, 2,3. 
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Further we describe in details the set MAX^ since it defines a sub-manifold in M with the 
maximum dimension 3, while MAX^ and MAX^ define sub-manifolds with dimension 1 and 2. 


Using new coordinates (58) in the case a 7 ^ 0 we get 

sn T cn r dn rfy (p) 


V = 


a\a\ce^k‘^{cn^ r — dn^ r sn 2 p)(l — k'^ sn^psn^ r) ’ 
fy{p) = —cn^pdnpsnp + (Ibae^/c^ E(p) — a'^p){dn^ p — k^ ct? psv? p). 

This gives us a description of the Maxwell set: 


MAX^niV^ = {z/ E Ar± I y{v) = o,r / 0,r G (-K,K),p G (0,K)} 


T± 


= {z^ G I fy{p) = 0,p G (0,K)}, 

where v = {t,p, E, a). 

The following obvious lemma will be useful for localization of roots of functions. 
Lemma 5. Let smooth functions f{u),g{u) satisfy on (0,rto) U M the conditions 


(59) 


u^O g(u 

Then f{u) >0 for u G (0, uq) ■ 


f(u) ^ 0, g(u) > 0, 

lim M=0. 


f(u) 

ff(u) 


> 0 , 


(60) 

(61) 


If functions / and g satisfy conditions (60), (61), then we say that g \s a. comparison function 
for / on the interval (0,uo)- 

Lemma 6. The function fy{p) > 0 forp G (0,K(A:)), k G (0,1). 

Proof We show that the function g{p) = A:^(l — A;^(l — cn^p)) is a comparison function for 
fy{p) for p G (0,K). 

The inequality fy{p) ^ 0 follows from the expansion fy{p) = ^a^k^p^ + o{p^). Notice that 
dip) > 0 for p G (0, K). Finally we get the equalities 

f fyip)\_ 4 a^sn^pcn^pdnp ^ fyjp) _ 4^3 .3^ 

^(p)/ 1 - A:2(l - cn^p) ’ g{p) 3 

So g{p) is a comparison function for /y(p); thus, it follows from Lemma that fy{p) > 0 for 
pG(0,K). □ 


Theorem 5. MAX^ = 
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Proof. Follows immediately from description of Maxwell set (59) and Lemma 


□ 


We proved that a general geodesic with A G has no Maxwell times tMAx(A) G (O, tsupr(A)) 
corresponding to the symmetries e*, i = 1,2, 3. So we conjecture that geodesics q{t) = Exp(A, t) 
with A G are optimal for t G (O, tgupr (^))- 


6 Spacelike normal extremal trajectories 

In the spacelike case H = K — + >0, we consider extremals on the level surface 

H = 1/2 and introduce new coordinates on this surface: 

hi = sinh 0 , /i 2 = ±cosh 0 , hs = c, /14 = a. 

Notice that we have the following symmetry of the Hamiltonian system: 

: (hi, /i 2 , hs, / 14 , xi, X 2 , y, z) ^ (hi, -h 2 , -hs, -hi, xi, -X 2 , -y, -z), (62) 

So we assume h 2 = cosh 0 > 0 without loss of generality in the sequel. 

In the variables {9, c, a, xi, X 2 , y, z) on the level surface {H = 1/2} the Hamiltonian system 
of Pontryagin maximum principle takes the following form in the normal case for spacelike 
curves: 


xi = — sinhh. 

(63) 

X 2 = cosh0. 

(64) 

X 2 sinh 0 + xi cosh 9 

y= 2 

(65) 

I 

• *^1 ' '^‘1 1 

2 : =-^^ cosh6<. 

( 66 ) 

9 = —c, 

(67) 

c = —a sinh 0, 

( 68 ) 

(5 = 0. 

(69) 


The vertical subsystem (67)"(69) has the energy integral: 

h2 (2 

E = — h 2 hi = -^— a cosh h. 


( 70 ) 


The family of all normal extremals is parametrized by points of the set 
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Figure 4; Phase portraits of the vertical subsystem for a = — 1, 1 (spacelike case) 


C = T*^n{H = l/2,h2 > 0} = {(hi,/l2,/l3,/l4) \ hl-hl = l,h2> 0} 

= {{9,c,a)eM.^}. 


The set C has the following decomposition into invariant sets of equations (67)“(69): 


C = Ci n Cj = 0, j, A = {6, c, a), 

Cl = {\ £ C \ a < 0,E > -a}, 

C 2 = {X £ C \ a > 0,E < -a}, 

C 3 = {X £ C \ a > 0,E > -a}, 

Ci = {X £ C \ a > 0, i? = —a, c 7 ^ 0, 0 7 ^ 0}, 

C 3 = {X £ C \ a ^ 0, E = —a, c = 0, 9 = 0}, 

Ce = {A G C I a = 0, E > 0, c ^ 0}, 

Cy = {A G C I a = 0, i? = 0, c = 0}. 

In order to parameterize extremals, we introduce coordinates (cp, E, a) on the set Ci 
in the following way. 
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Figure 5: Trajectories of the vertical subsystem for a = —1, 1 (spacelike case) 


In the domain Ci: 
k = 


E cy. 


e (0,1), 


E — a 
c = \/2{E + a) sn(as(/j), 

\/E — a y/E + acn(8e(p) dn(8e(p) 


sinh0 = 

In the domain C 2 ' 
k = 


a 


2 a 


G (0,1), 


c = sgn 9 
sinh 9 = sgn 9 
In the domain C 3 : 


a — E 

Y^2(—a — E) sn(as(p) 


cn(ae(p) ’ 

\/—a — E \/a — E dn(8e(p) 
acvL^{?e^p) 


k = ^ k 2 = k^ & (- 00 ,1), 


c = sgn o 
sinh0 = sgnc 


E a 

2asdn(8e(/9) 


cn(8e(/9) 
2 sn(ae(p) 
cn 2 (ae(p) ’ 


Remark 2. If k'^ < 0, then k = i 


k = 


ae = 


(p £ 


E — a 


cosh 9 = 1 — 


(E + a) ciF(as(p) 


a 


as = 


a — E 


K K\ 
ae m/ 


cosh 0 = 1 — 


a -\- E 
a cn(ae<p) 


m = 


F G - 


E + a 
2 ’ 
K K 


0S 96 


cosh0 = 


1 + sn^(as(p) 
cn2(as(p) 


\/l-fc2’ Vl-/c2 


-fc 2 


G (0,1). Using formulas from we 
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have 


k = 

sn(ae(/?, k) = 


a — E 


l-k^ = 


2a ’ 

sn{^/a^p, k) 


a + E 


2a 


dn(\/a(/?, k) 


dn(8ev9, k) = 
cn(8ev9, k) = 


dn(-y/a(/j, k) ’ 
cn(^(/?, fc) 
d.T^{^/aip^ k) 


(71) 

(72) 


In the domain C^: 

c = sgn f 
sinh 9 = sgn ( 


I_ 

Q,e2a8</j _ 1 ’ 

4me®‘^(ae^®‘^ + 1 ) 

- 1)2 ■ 


ae = -v/a) G M, 7 I 1 , 

ae2®¥’(ae2*‘^ + 6) + 1 


cos 0 = 


ae 


2 ae(/? _ 


Immediate differentiation shows that in the coordinates {ip, E, a) the subsystem for the 


costate variables (67|-(69) takes the following form: 


ip = I, E = 0, d = 0, 


(73) 


so that it has solutions 


ip{t) = (ft = ifo + t, E = const, a = const. 


(74) 


6.1 Exponential mapping for spacelike normal extremals 


Denote ^po = seipo,'tpt = se<pf From the definition of the variables c and 0 we obtain the 
following parametrization of extremal trajectories. 

If A G Cl, then 


Xl{t) 

X2{t) 

y{t) 

z{t) 


\f2\J a + E 


a 


- snV’o), 


y/a + E 

v/ 2 a 2 


- t. 


2 ae( cn ipt dn ipt - cn-ipodmpo + (E(' 0 t) - E(V'o)) (sn ipt + sn 


+ al^snipt + sn tj , 

(X fill ^ 1 /" 

28e(^(a + E)(cnV’tdnV’t(snV’t -sni/'o) 

— 2{cn'ijjt dn'ijjt — cnipo dn'i/ig) sn?/)o) — {E + 3(a + E) sn^ V’o)(E(V'i) 
+ a(^a — E — 3(a + E) sr 


(75) 

(76) 


(77) 


(78) 


Projections of extremals to the plane {xi,X 2 ) with \a\ = 1 and pJo = 0 are shown in Fig. left. 
This case has no upper bounds for time: tsupr(A) = + 00 . 
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Figure 6 ; Projections of goedesics to the plane {xi,X 2 ) for A G Ci, A G (^2 and A G (Pa 
(spacelike case) 


If A G 6 * 2 , then 

y/2y/—a — E sgn 6 


xi{t) = -- 
X2{t) = - 


a 


(^scipt — sc 


2as(E(fo) - E('!/;o) + dnV’oscV’o - dnV’t sc V't) + Et 


a 


(79) 

(80) 


yit) = - 


y/—a — E sgn 9 


V2a^ 


E sc '!/^o + sc V't) t + 236 (^ (sc -il^o + sc Tpt) (E(V't) - E('!/;o)) 

+ (dnV’t-dnV'o)(l-scV'osc'^t)) ), ( 81 ) 


u\ (®2(t))^ , 2as(3(a + E)sfo'!/)o-EcfoV>o), . . 

z{t) = -+ ————-o ^ E{^Pt) - 


+ 


3q:^ cn^ ipo 

2 as| cn'0d q;( cfo'04 dn^^ 0(1 ~ 3cn^ ^/^o) 


30^38 cn^ tpt cfo 00 

+ 3 cfo 00 dn ipt) + cfo 0t dn 00 (1 — 4 cn^ 0o) +3 cn^ 00 dn iptj E ) sn 0 o 
— cn 00 dn 0t I a cn^ 0o + 3a sfo 0o cn^ 0t + ( 3 cfo 0t + cn^ 0o (l 


- 4cfo0t) )E sn0t + 


(a + E) (cfo 00 (a — 4E) + 3E) 


(82) 


3a^ cn^ 00 

Projections of extremals to the plane (xi,X 2 ) with |a| = 1 and 0o = 0 are shown in Fig. 
center. Notice that since ipt G (—K,K) we have the upper bound for the time parameter: 

K{k) - 00 


t ^ tsupr(A) — 


36 


(83) 
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If A G Cs, then 


Xl{t) 

X2it) 

y{t) 

z{t) 


2 sgn c(cn ^po dn — cn tpt dn tpo) 

ae cn'f/jj cn V’o (1 — ^ 2 ) 


3e(l - k 2 ) 
sgnc 


dn ^lJt sc ipt - dn V’o sc V'o - E(V't) + 
2 


+ 1, 


■(^{dc'tjjt + dcV’o) (E(V^t) - E(V'o)) 


ae(l - k 2 ) V®(1 “ ^ 2 ) 

+ (dc'i/’o dc V't + ^2) (sn ?/)o -sn-i/’t)) - (^dcV’o + dcV’t^t^ , 

2{E{ipt) -E(V;o))(6-6/c2 + cn^V>o(l + 7fc2)) 2(3dn^ V^p + cn^ ipp) 


+ 


3as^ cn^ ipo{—l + ^ 2 )^ 3®^ cn^ V’o(—1 + ^ 2 )^ 

cn^ V’t dn 'i/>o (2 + cn^ V'o — 2 A ;2 + 7 cn^ V’ 0 ^ 2 ) 


386^ cn^ cn^ V’o(~l + ^ 2 )^ 

X sn?/)o — 6 cn^ -ipo cndn -ipo (/c 2 — l) sn'tpt + 2 cn^ 'tpQ dn %l)t {k 2 — l) sn 

( 

- cnV’pcn^ ^|;tdn^|;t{6 - 6k2 + cn^ V'o(l + 7^2)) snV’t^ H- ^ -• 


(84) 

(85) 


(86) 


(87) 


Projections of extremals to the plane (xi, X 2 ) with \a\ = 1 and ■00 = 0 are shown in Fig.|^ right. 
Notice that since ^pt G (— K, K) we have an upper bound for time parameter: tsupr(A) = , 

If A G (74, then 


/ gbt gbo \ 

^ l+ae^^ —l + ae^’^J ’ 

, , ^2e’^o(-l + ae’^*+4^o)(2 + ast + - 2)) 

y{t) = -sgne^— - / 1 ^ 2^a\r 1 ^ 2 ^a - 

ae(—1 + -g 

x^{t) 4 / —l — 9ae‘^^°{—2 + ae‘^^°) 

6 3 Y (—1 + ae^’^o)^ Q;8e(—1 + ae^’^o)^ 

12(1+ ae2V'o(-l+ 2e4/’‘-^o)) 8 

Q:as(—1 + ae^^°){ — l + ae^'^^Y Q:as(—1 + 

3(^1 + ae^Po (6 + - ae2Po(-l + 4e^^-'l’°))^ 

Q:8e( —1 + ae'^'^°Y{—l + ae^^*) 


( 88 ) 

(89) 

(90) 


(91) 


Projections of extremals to the plane (xi, X 2 ) with \a\ = 1 are shown on the left (for 0o<O) and 
in the center (for '0o<O) in Fig. 

Notice that here we have the condition Y 1- Denote tp =- 1 =~—■ If tp < 0, 

2Ya - (/?p 

then there is no upper bound for time parameter: tsupr(A) = +00 (see Fig. left). If tp > 0, 
then tsupr(A) = tp (see Fig. center). 
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Figure 7: Projections of goedesics to the plane {xi,X 2 ) for A G C 4 (2 types) and A G Ce 
(spacelike case) 

In other cases the Hamiltonian system is easy to integrate. If A G 6 * 5 , then 

xi{t) = 0, X 2 {t) = t, y = 0, z = —. (92) 

0 

Projection of the extremal to the plane (xi,X 2 ) is the vertical straight line. This case has no 
upper bounds for the time: tsupr = + 00 . 

If A G Ce, then 


Xl{t) 

X2{t) 

y{t) 

z{t) 


0 = 0 , c = const / 0 , 6 

cosh(ct — Oq) — cosh 00 
c 

sinh(ct — 0 o) + sinh 0 o 


ct — sinh(ct) 

2? ’ 

4cosh(^-30o) sinh^f- 


6c3 


00 - ct, 


3 cosh 00 (ct 


— sinh(ct)) 


(93) 

(94) 

(95) 

(96) 


Projections of extremals to the plane (xi,X 2 ) with |q;| = 1 and ^jJo = 0 (hyperbolas) are shown 
in Fig. right. This case has no upper bounds for time: tsupr = + 00 . 

If A G C 7 , then 


a = 0 , c = 0 , 0 = const, sinh 0 = sq, cosh 0 = co, 


xi(t) = -Sot, 


X2{t) = Cot, 


y{t) = 0, 
z{t) = 


(97) 

(98) 

(99) 
( 100 ) 
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Projections of extremals to the plane (xi,X 2 ) are straight lines. This case also has no upper 
bounds for time: tgupr = +oo. 

So we obtained parameterization of the exponential mapping 

Exp : iV —)• M = Exp(A,t) = q{t), (101) 

N = {(A,t) € C X M_|_ I t G (O, tsupr(A))}, (102) 

in the spacelike case. It maps a pair (A,t) to the point of corresponding extremal q{t). 

Eurther we investigate discrete symmetries of the exponential mapping and on this basis 
find estimates for the cut time on extremals. 


6.2 Discrete symmetries of exponential mapping 


Subsystem (67)-(69) for the costate variables of the normal Hamiltonian system has sym¬ 
metries which preserve the direction field of this system. 




Figure 8: Action of symmetries on trajectories of vertical subsystem for a = —1,1 
(spacelike case) 


Let us describe the action of symmetries e* on the set of trajectories of the vertical subsystem 
with preservation of time direction. Denote a smooth curve 

7 = I t G [0,r]| C C. 
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Define action of symmetries on these curves (see. Fig. [^: 


: 7h^ 71 = ^{ 6 ^(t), {t), a^) \ t G [0,r]| = | ((9(T - t), -c(r - t),a) |, 

: 7 ^ 72 = |c^(t), a^) | t G [0, r]| = | ( - e(T - t),c{T - t),a)'^, 

: 7 ^ 73 = | (6'^(t), c^(i), a^) I i G [ 0 , F]| = | ( - e{t), -c(t),a)|. 


It is obvious that if 7 is a solution to the vertical subsystem (67)-(69), then 77 f = 1,2,3, is a 
solution as well. 

Consider the group of symmetries G = {Id, o e^} = 7? x Z^. The symmetry 

preserves direction of time, while and change it. 

Continue the action of symmetries from the vertical subsystem to solutions of the Hamilto¬ 
nian system of PMP 


6{t) = —c{t), c{t) = —asmh6{t), a = 0, 
q{t) = - cosh 9{t)XI [q{t)) + sinh 6{t)X 2 {q{t)), 


in the following way: 

e* : |(6»(t),c(f),a,g(f)) | t G [0,r]| ^{e\t),d{t),a\q\t))\t G [0,T]|. (103) 

Let q{t) = X 2 {t),y{t), z{t)), t G [0, T], be a geodesic and 

q\t) = {x\{t),xl{t),y\t),z\t)), t G [0,r], i = 1,2,3, 
be its image under the action of a symmetry e*. 

Lemma 7. The symmetries e* map trajectories in the plane (xi,X 2 ) in the following way: 

x{(t) = Xl(T) - Xl(r - t), xl{t) = X2{T) - X2{T - t), 

xl{t) = Xl{T - t) - Xl(r), X^it) = X2{T) - X2{T - t), 

xl{t) = -Xl{t), xl{t) = X 2 {t). 


Proof. Using direct integration, for instance, for we verify that 

®i(^) = f { — si'nh9{T — s))ds = — f [ —sinh9{r))dr = xi{T) — xi(T — t), 

Jo Jr-t 

® 2 (^) = [ cosh 9(T — s)ds = — f cosh 9{r)dr = X 2 {T) — X 2 {T — t). 

Jo Jr-t 

□ 
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Lemma 8. The symmetries e* maps endpoints of geodesics q = [xi,X 2 ,y, z'j to endpoints of 
5 * = (x^, X 2 , y*, 2*) in the following way: 

xl(T)=xi(r), x^(r) = X2(r), y\T) = -y{T), z\T) = z{T) - xi{T)y{T), 

x?(r) =-xi(r), xi(r) = X2(r), y\T)=y{T), z\T) = z{T)-xi{T)y{T), 

x?(r) =-xi(r), xi(r) = X2(r), y\T) = -y{T), z\T) = z{T). 

Proof Lemmagives us expressions for x\{T) and x^iT). The expressions for other variables 
are obtained by integration. For example, for we have 

(D = £ 4(t}sinheHt)+^4{t)cosbsHt)^^^ = ^ 






10 


2 G ^ Jo 

{ x 2 { t ))\ /■'^(xi(r))'-2xi(r)xi(r-t) + (xi(r-t))\ 

- 2 -^ io - 2 - -tjdt 

(MT)f , (ii(r))"i2{T) /-r (x2{T)f 

H-;;-a;i(r) / xi(s)x 2 (s)as + 2 ;(r)- 


6 


6 


- MT) {y{T) + ^ ^ ^ _ 3;^(T)y(r). 


□ 

We define the action of e® in the preimage of exponential mapping N by restricting action 
to the initial point of trajectory of the vertical subsystem: 


e^:C^C, e\e,c,a) = {e\c\a^), 
{6^, c^, a^) = (0, —c, —a), 

(6»2,c^a2) = (-6»,c,a), 

( 0 ^,c^,q:^) = (— 0 ,—c, —a), 


(104) 

(105) 

(106) 
(107) 


in the following way: 


e®(A,t) = (e® o e*^’'(A),t), i = l,2, 
e^(A,t) = (e^(A),t), 
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where — asinh0^ G Vec(C') is the vertical part of the Hamiltonian vector field. 

We define the action of e* in the image M of the exponential mapping M by restricting the 


action to endpoints of geodesics (see Lemma [^: 

e" : MM, e\q) = e\xi,X 2 ,y,z) = = {x\,xl,y\z"), (108) 

{x\,xl,y^,z^) = {xi, X 2 , -y, z-xiy), (109) 

{xl,xl,y^,z‘^) = {-xi, X 2 , y, z - xiy), (110) 

{xi,xl,y^,z^) = {-xi, X 2 , -y, z). (Ill) 


Since the actions of e* in the domain N and the image M of the exponential map are induced 


by the actions of symmetries (103) on trajectories of Hamiltonian system, we have the following 
result. 


Proposition 2. The mappings = 1,2,3, are symmetries of the exponential mapping, i.e., 
(e* o Exp)(0, c,a,t) = (Exp oe^){6, c,a,t), {6,c,a,t) £ N. 


6.3 Maxwell points of spacelike normal extremals 

In this section we compute the Maxwell points corresponding to the symmetries e*. On this 
basis we derive estimates for the cut time along extremals. 

We define the Maxwell sets in the preimage of Exp corresponding to the symmetries e*: 

MAX* = {(A,f) G I AV A, Exp(A\t) =Exp(A,t)}, A = (0,c,a), A* = e*(A). 

It follows from Proposition that the equality Exp(A*,t) = Exp(A,t) is equivalent to 
£*(9(^)) = Qi't)- Therefore we get a description of fixed points of the symmetries e* in the 
preimage of the exponential mapping. 

Lemma 9. 


1. e^{q) = q y = 0. 

2. e^{q) = q xi = 0. 

3. e^{q) = q xf + = 0. 


Proof. Follows from definition (108)-(111) of action of the symmetries in M. 


□ 
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Now we compute fixed points of symmetries in the preimage of the exponential mapping, 
which is necessary for description of Maxwell sets. Introduce the following coordinates in the 
set Ni = G CiXR+ \t G (O, tsupr(A))}: 


r = 


V’t + V’o 


P = 


il^t - -00 ast 


( 112 ) 


1. = A 


2 ’ " 2 2 
The parameter r corresponds to the middle point of a trajectory. 

Lemma 10. 

/ 

snr = 0 , z/AgC'iUC 2 , 
impossible, z/ A G Cs U (74 U Ce, 
c = 0, if X G U Ct- 

cnr = 0, if X G Cl, 
impossible, if X G C 2 U C 4 , 
snr = 0, if X G C 3 , 

e = o, tfXGuUCi- 

impossible, if X G (Ui^=i ^i) U Cq, 

9 = c = 0 ,ifXGC 5 UC 7 . 

Proof. Follows from definition of action of the symmetries in the domain of the exponential 

□ 


2. A^ = A 


3. A^ = A 


mapping (104)-(107) and definition (112) of the coordinate r. 


The equations e*(g) = q define sub-manifolds in with dimension from 3 to 2, the expo¬ 
nential mapping transforms the corresponding Maxwell sets into these sub-manifolds: 

Exp(MAX*) C {g G I e\q) = q} i = 1,2, 3. 

Further we describe in detail the sets MAX^, MAX^ since they define sub-manifolds in M 
with the maximum dimension 3, while MAX^ defines a sub-manifold with dimension 2. We 


investigate roots of the equations xi = 0, y = 0 along a geodesic. Using coordinates (112) for 
the case A G (7i we get 

2^2{a P E) cn r dn T snp 


Xl = 


y = 


q:( 1 — sn^psn^ r) 
2 y/2{a + E) sn r/i (p) 


a^ae(l — kf sn^psn^ r) ’ 
flip) = apcnpdvip + as^(2 cnpdnpE(p) — (1 -|- k^) snp -|- 2k^ sn^p). 
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This gives us description of the Maxwell sets in N: 

MAX^ nA^i = {i^ € Ni \ y = 0,snr / 0} = {i^ G | fi{p) = 0}, (113) 

MAX^ nA*! = G Ai I xi = 0, cnr / 0} = {j/ G Ai I snp = 0}, (114) 

where ly = {t,p, E, a). 

Let p\,i £ N, denote a positive root of equation fi{p) = 0, s.t. p\ < pf < ■ ■ ■ < p\ < ■ ■ ■ ■ 
We have iMAXi("^) ~ -Pi) ^max2 ("^) = 2K. There is a conjecture (this conjecture is supported 
numerically) that tMAX^("^) < ^MAX^("^)- This means that each geodesic from Ci meets the set 
MAX^ hr St. 

If A G (72, then 

2 sgn —a — Ednrcnpsnp 

7 ^ ^ j 2 A ’ 

a(cn^ r — sn^pdn r) 

_ 2sgn0\/2\/ —a — E CUT surf 2 {p) 

a^8e(cn2 T — sn^pdn^ r) ’ 

f 2 {p) = —E{p) dnp — ae^ (2 dnpE(p) — A:^ cnpsnp) (115) 

We have 

MAX^nAa = {z/ G Aa I y = 0,snT / 0,r G (-K,K),p G (0,K)} 

C {z/ G Aa I/2(p) = 0,p G (0,K)}, (116) 

MAX^nAa = {z/ G Aa I XI = 0,T G (-K,K),pG (0,K)} = 0, (117) 


where v = {t,p, E, a). 

Lemma 11. The function f 2 {p) > 0 for p G (0,K(A:)), k G (0, 1). 


Proof. We show that g{p) = dnp is a comparison function for /2(p) for p G (0, K). 

The inequality f 2 {p) ^ 0 follows from the expansion f 2 {p) = \ak‘^p^ + o{p^). Notice that 
g{p) > 0 for p G (0, K). Finally we get the equalities 


( f2{p) \ 
V 9{P) ) 


2 2 2 

a sn pen p 

se^ dn^ p ’ 


^ = latV + c(p=>). 


So g{p) is a comparison function for /2(p); thus, it follows from Lemmathat f 2 {p) > 0 for 
pG(0,K). □ 


Proposition 3. MAX^ PlAa = 0 . 
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Proof. Follows immediately from Lemma 11 and (116). 


If A G C 3 , then 

4sgncsnpsnr 
as(cn2 r — dn^ r sn^ p) ’ 

4 sgn c cn r dn r/s (p) 

^ ae2(l —/c2)^(cn2 T — dn^ r sn^p) ’ 

f^{p) = cnpdnp((l — k 2 )p — 2 E(p)) + (1 + ^2) snp — 2/c2 sn^p. 


□ 


(118) 


We have 


MAX^nlVa = {i^ G iVs I 2/ = 0,r G (-K,K),p G (0,K)} 

= e I/3(p) = 0,p G (0,K)}, (119) 

MAX^nlVa = {i^ G iVs I XI = 0,snT 7 ^0,r G (-K,K),p G (0,K)} = 0, 


where ly = {t,p, E, a). 

Lemma 12. The function fz{p) > 0 forp G (0,K). 


Proof. We show that g{p) = cnpdnp is a comparison function for fsip) for p G (0, K). 


(y.^ o o 

The inequality f^ip) ^ 0 follows from the expansion fslp) = — jp + oip ). Notice that 
g{p) > 0 for p G (0, K). Finally we get the equalities 

(h{p)\ _ {l-k2fsn?P f2{p) _ «\3 I 3^ 

\ g{p) } cn^pdn^p ’ g{p) 

So g{p) is a comparison function for fsip)', thus, it follows from Lemmathat f 3 {p) > 0 for 
p G (0, K). □ 


Remark 3. Notice that if k"^ = k 2 < 0, then Lemma 12 is valid. We have the following 
formulas: 


xi = 


y = 


2y/2\Ja + E sgncdnr snrdnpsnp 
a(cn2 r — dn^ r sn^p) ’ 

2\/2 cn TyJ a + E sgn c/4(p) 

am(cn2 r — dn^ r sn^ p) ’ 


ae = m = ya 


/4(p) = 2cnpE(p) — pcnp — dnpsnp, k = k = 


a — E 
2a 


One can easily prove that g{p) = cnp is a comparison function for fi{p) for p G (0, K). 

Proposition 4. MAX^ ON^ = 0. 
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Proof. Follows immediately from Lemma 12 and (119). 


□ 


We found the following upper bound for the cut time along general extremal with A G 




Theorem 6. 


A G Cl 

fcut(A) < 2K(fc). 

( 120 ) 

A G C 2 U C 3 = 


( 121 ) 


There is a conjecture that this bound is sharp. 


7 Conclusion 

In this paper we consider extremals in the Engel group with a sub-Lorentzian metric. For 
this problem 

• we prove that timelike normal extremals are locally maximizing; 

• we calculate the timelike, spacelike and lightlike normal extremal trajectories; 

• we describe the symmetries of the exponential map and the corresponding Maxwell points 
of timelike, spacelike normal extremals; 

• we obtain an upper estimate for the cut time on extremal trajectories. 

In the future we are intending to investigate the optimality of extremal trajectories, relying 
on the method of nonlinear approximation we shall also apply our results to consider general 
non-linear systems with 2-dimensional control in a 4-dimensional space with a sub-Lorentzian 
metric. 


References 

[1] A.A. Agrachev, El-H. Chakir El-A.,and J. P. Gauthier, Sub-Riemannian metrics on 
Proc. Conf. Canad. Math. Soc. 25(1998). 

[2] A. Agrachev, D. Barilari, and U. Boscain, On the Hausdorff volume in sub-Riemannian 
geometry. Calculus of Variations and Partial Differential Equations 10.1007/s00526-011- 
0414-y (2011), 1-34. 


35 





[3] A. Agrachev, U. Boscain, J.-P. Gauthier, and F. Rossi, The intrinsic hypoelliptic Laplacian 
and its heat kernel on unimodular Lie groups, J. Fund. Anal. 256(2009), No. 8, 2621-2655. 

[4] A. A. Agrachev, G. Gharlot, J. P. A. Gauthier, and V. M. Zakalyukin, On sub-Riemannian 
caustics and wave fronts for contact distributions in the three-space, J. Dynam. Control 
Systems 6 (2000), No. 3, 365-395. 

[5] A. Ardentov, Yu. Sachkov, Extremal trajectories in a nilponent sub-Riemannian problem 
on the Engel group, Matematicheskii Sbornik 202(2011), No. 11, 31-54. 

[6] A. Ardentov, Yu. Sachkov, Gonjugate points in nilpotent sub-Riemannian problem on the 
Engel group, J. Mathematical Sciences 195(2013), No. 3, 369-390. 

[7] A. Ardentov, Yu. Sachkov, Gut time in sub-Riemannian problem on Engel group, ESAIM: 
Control, Optimisation and Calculus of Variations, accepted (2015) 

[8] N. I. Akhiezer. Elements of the theory of elliptic functions. Providence, RI: American 
Mathematical Society. 1990. 

[9] Gartan, A, Sur quelques quadratures clout I’elernent diflerentiel contient des fonctions 
arbitraires. Bull. Soc. Math. France 29(1901), 118-130. 

[10] Q. Gai, T. Huang, Yu. Sachkov, X. Yang, in the Engel group with a sub-Lorentzian metric, 
submitted. 

[11] El-H. Alaoui, J-P. Gauthier, and I. Kupka, Small sub-Riemannian balls in M^, J. Dynam. 
Control Sys. 2(1996), No. 3, 359-421. 

[12] R. Beals, B. Gaveau, and P.G. Greiner, Hamilton-Jacobi theory and the Heat Kernal on 
Heisenberg groups, J. Math. Pures Appl. 79,7(2000) 633-689. 

[13] J.K. Beem, P.E. Ehrlich, and K.l. Easley, Global Lorentzian geometry, Marcel 
Detter(1996). 

[14] D. G. Ghang, I. Markina, and A. Vasiliev, Sub-Lorentzian geometry on anti-De Sitter 
space, J. Math. Pures Appl. 90(2008), No. 1, 82-110. 

[15] M. Golubitsky and V. Guillemin, Stable mappings and their singularities, Spinger-Verlag, 
New Yor^(1973). 


36 



[16] M. Grochowski, Differential properties of the sub-Riemannian distance function, Bull. Pol¬ 
ish. Acad. Sci. 50(2002), No. 1. 

[17] M. Grochowski, Geodesics in the sub-Lorentzian geometry. Bull. Polish. Acad. Sci. 
50(2002), No. 2, 162-178. 

[18] M. Grochowski, Normal forms of germs of Gontact sub-Lorentzian structures on M^, Dif¬ 
ferentiability of the sub-Lorentzian distance function, J. Dynam. Control Sys. 9(2003), No. 
4, 531-547. 

[19] M. Grochowski, Reachable sets for the Heisenberg sub-Lorentzian structure on M^, an 
estimate for the distance function, J. Dynam. Control Sys. 12(2006), No. 2, 145-160. 

[20] M. Grochowski, On the Heisenberg sub- Lorentzian Metric on R^, Ceometric Singularity 
Theory, Banach Center Publications, 65(2004), 57-65. 

[21] M. Gromov, Garnot-Garatheodory spaces seen from within, Progr. Math. 144, Birkhauser, 
Boston (1996), 79-323. 

[22] T. Huang, X. Yang, Geodesics in the Heisenberg Group Hn with a Lorentzian metric, J. 
Dynam. Control Sys. 18(2012), No. 1, 21-40. 

[23] F. Monroy-Perez, and A. Anzaldo-Meneses, Optimal control on the Heisenberg group, J. 
Dynam. Control Sys. 5(1996), No.4, 473-499. 

[24] B. O’Neill, Semi-Riemannian Geometry: with Applications to Relativity, Pure and Applied 
Mathematics, vol. 103, Academic Press, Inc., New York, 1983. 

[25] A. Korolko, and 1. Markina, Non-holonomic Lorentzian geometry on some H-type groups, 
J. Ceom. Anal. 19(2009), 864-889. 

[26] R. Montgomery, Singular extremals on Lie groups. Math. Control, Signals and Systems 
vol. 7(3) (1994), 217-234. 

[27] R. Montgomery, A tour of sub-Riemannian geometries, their geodesics and applications. 
Math. Surveys and Monographs 91, American Math. Soc., Providence, 2002. 

[28] P. Piccione, and D.V. Tausk, Variational aspects of the geodesic problem in sub- 
Riemannian geometry, J. Ceometry and Physies. 39(2001), 183-206. 


37 



[29] Yu. L. Sachkov, Conjugate and cut time in the sub-Riemannian problem on the group of 
motions of a plane, ESAIM Control Optim. Calc. Var.16 (2010), 1018-1039. 

[30] R. Strichartz, Sub-Riemannian geometry, J. DijJ. Ceom., 24(1986), 221-263. 


38 



